Annular Supersonic Swirling Flows with Heterogeneous Condensation 
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In recent years, separating and extracting technologies of condensate gas have been developed by combining a 


swirl flow with non-equilibrium condensation phenomena of condensate gas generated in a supersonic flow. The 


technology can reduce the size of the device and does not use chemicals. However, there are many unresolved 


problems for performance of the separation, extraction and operating principle. Therefore it is necessary to re- 


search further in order to improve the performance of the equipment. In the present study, the numerical study 


was carried out to clarify the effect of the heterogeneous condensation on the characteristics of the swirling flow 


field in a supersonic annular nozzle, and the differences between homogeneous condensation and heterogeneous 


condensation in the flow field. As the results, it is found that the condensation flow with a swirl affects the posi- 


tion of sonic line, the generating position of condensate and the radial distribution ratio of liquid phase. 
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Introduction 


Condensation phenomenon of moist air or steam in a 
high speed flow field is, in general, induced through the 
non-equilibrium condensation process [1]. In this process, 
condensation nuclei are generated by collision and ag- 
gregation of vapour molecules. The condensation of the 
vapour, so called homogeneous condensation, takes place 
on the nuclei. On the other hand, in heterogeneous con- 
densation [2], the condensation of the vapour takes place 
on foreign nuclei : smoke and vapor from fires and vari- 
ous industries, dust from land surfaces, salt from ocean 
and particulate products from chemical reaction. Further, 
effects of foreign nuclei on the high speed flow field 
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have not been clarified significantly. 

Recently, we are faced with solution of environmental 
problems and the sophistication of the energy structure 
on a global scale. Natural gas of lower CO, emissions is 
important as an alternative energy of petroleum and coal 
[3]. The separation and extraction of condensable gas and 
impurities contained in the natural gas are important for 
natural gas processing. Separating and extracting tech- 
nologies of condensable gas have been developed by 
using a swirling flow in a supersonic annular nozzle [1, 4, 
5]. The technology can reduce the size of the device and 
it has an advantage that does not use chemicals. However, 
there are many unresolved problems for performance of 
the separation, extraction and operating principle. There. 
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Nomenclature 


Q Conservation mass term Greek symbols 

EEG Numa Ms E ae 
E,, F,,G, Source term of viscosity o* Displacement thickness, (m) 
I Source term of turbulence ġ Relative humidity, (%) 

S Source term of condensation p Density, (kg/m*) 

D Diameter, (m) v Kinematic viscosity, (m’/s) 

g Condensate mass fraction oO Specific dissipation rate, (1/s) 
I Nucleation rate per unit time and volume, (1/(s:m*)) Subscripts 

k Specific turbulence kinetic energy, (J/m’) 0 Stagnation point 

L Latent heat, (J/kg) v Vapor 

n Number density, (1/kg) a Dry air 

Nhet Number of embryos per unit mass, (kg”) e Experiment 

p Pressure, (Pa) het Heterogeneous 

Re Reynolds number hom Homogeneous 

r Radius, (m) in Inner 

Ry Radius of a solid particle, (m) l Liquid 

Swain Swirl number out Outer 

T Temperature, (K) ee) Infinite flat plane 

u, v, W Velocity components, (m/s) * Dimensionless variable 
Xyz Cartesian coordinates, (m) 


fore it is necessary to research further in order to improve 
the performance of the equipment In the previous re- 
search, by using a homogeneous condensation pheno- 
menon of moist air occurred in the supersonic flow in the 
annular nozzle with a swirl, the effect of the nozzle inlet 
shape on spatial distribution of condensate was investi- 
gated numerically, and the possibilities of the condensa- 
ble gas separation and reduction of device size was 
shown qualitatively [6]. However, the effect of heteroge- 
neous condensation on the supersonic annular flow with 
a swirl has not been taken into account in the study. 

In the present study, the numerical investigation is 
carried out to clarify the effect of the heterogeneous 
condensation on the characteristics of the swirling flow 
field in a supersonic annular nozzle and the differences 
between homogeneous condensation and heterogeneous 
condensation in the flow field. As the results, it is found 
that the condensation flow with a swirl affects the posi- 
tion of sonic line, the generating position of condensate 
and the radial distribution ratio of liquid phase. 


Experimental Apparatus and Method 


In order to validate the present computation, an expe- 
rimental work was conducted for homogeneous conden- 
sation phenomenon. Figure 1 shows a schematic diagram 
of experimental apparatus. A supersonic indraft wind 
tunnel where the moist air at atmospheric pressure is 


drawn into a vacuum tank, was used in the present expe- 
riment. 

Figure 2 shows details of test section. An annular noz- 
zle is composed of an inner body and an outer nozzle. 
The outer nozzle is coaxial with the inner body. Throat 
diameter of outer nozzle in case of no inner body is De = 
30 mm and design Mach number at exit in the annular 
nozzle is 3.0. The pressure measurements were con- 
ducted by a pressure transducer (Kulite XT-190) mounted 
on the pressure holes (diameter : 1.0 mm) along the outer 
wall. The supersonic steady flow lasting about 30 se- 
conds was obtained in the test section. The stagnation pre- 
ssure po and the stagnation temperature To of moist air in 
the reservoir were set at 101.2 kPa and 293 K, respective- 
ly. The values of relative humidity ø are 25% and 70%. 
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Fig. 1 Experimental apparatus 
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Fig. 2. Supersonic nozzle with inner body (D. = 30 mm) 


Assumptions in the present calculation are as follows ; 
both velocity slip and temperature difference do not exist 
among condensate droplets, solid particles and inert gas 
mixture, and the effect of the condensate droplets and 
solid particles on pressure is neglected. All particles are 
assumed to have a smooth and spherical shape, and all 
condensation nuclei are assumed to be chemically inert 
and insoluble in water vapor. 

The governing equations used in the present study are 
three dimensional compressible Navier-Stokes equations 
combined with equations of continuity, energy, turbulent 
kinetic energy, specific dissipation rate, conservation of 
mass of the liquid phase and conservation of the number 
density of droplets with homogeneous and heterogeneous 
nucleations. To estimate the eddy viscosity, k-@ model 
was employed in computations [7]. Non-dimensional 
governing equations are written in the three-dimensional 
Cartesian coordinate system (x,y,z) as follows : 

ðQ OE OF ðG 
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In Eq.(1), Q is conservative vector, E and F, G are in- 
viscid flux vector and E,and F,, G, are viscous flux vec- 
tors. J and S are the source terms corresponding to turbu- 
lence and condensation, respectively [8]. 

Total flow rate of liquid phase g [8] is written as 


E = Shom t Shet (2) 

where 2hom and Zne are condensate mass fractions 

generated by homogeneous and heterogeneous condensa- 

tions, respectively. Both of the condensate mass fractions 

g are expressed as a rate equation, based on the following 
equation [9]. 


Oe L 
dt = 3 Pl r(t) Pan (t) O) 
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The density of liquid phase p(T) is also a function of 
temperature [10] and the density of gas mixture p,, is 
calculated by the sum of density of vapor p, and dry air 


Pa r and r, are averaged and critical droplet radius, respec- 
tively. Z is nucleation rate per unit time and volume. 

Heterogeneous nucleation process consists of four 
nucleation stages [8]; I. Generation of embryo, II. 
Growth of embryo, III. Formation of liquid film and IV. 
Growth of liquid film. In the first stage, embryos are 
generated on the surface of a solid particle. In the second 
stage, embryos on the surface of a solid particle grow up 
and new embryos are generated on the particle. In the 
third stage, a liquid film is formed on the whole surface 
of a solid particle and the solid particle covered by the 
liquid film can be assumed a nucleus. In the fourth stage, 
the liquid film grows up. 

There are two models for heterogeneous nucleation 
process [8]. For model 1, four nucleation stages from I to 
IV are considered for nucleation process. For model 2, 
only the fourth stage (IV) of nucleation process is used. 
This model is assumed that nucleation stages from I to III 
proceed in an infinitesimal time. In the present study, 
model 2 is used as a heterogeneous condensation model 
because the difference between results obtained by both 
models is very small. Equations for heterogeneous con- 
densation in the present simulation were given based on 
Refs.[8,11,12]. 

Equations were discretized in space using a cell- 
centered finite volume formulation with a quadrilateral 
structured cell system and in time using the Euler impli- 
cit method [13]. Furthermore, a third-order accurate 
MUSCL TVD schemes based on Roe’s approximate 
Riemann solver [14] were applied to the inviscid fluxes 
and the viscous fluxes were evaluated by second-order 
accurate central differences method. Generation terms 
associated with turbulence and condensation were eva- 
luated by first-order accurate central differences method. 
Unfactored implicit equations derived with no approx- 
imate factorization were solved by a point Gauss-Seidel 
relaxation method [15]. 


Computational Conditions 


Figure 3 shows nozzle geometry used in the present 
study. An annular nozzle used in this calculation is com- 
posed of an inner body and an outer nozzle. The outer 
nozzle is coaxial with the inner body. Throat diameter of 
the outer nozzle in case without the inner body is D = 10 
mm which is the same as one at the nozzle exit. The noz- 
zle opening angle a is 4°. Design Mach number for the 
nozzle is 2.24. 

Figures 4 (a), (b) and (c) show a computational do- 
main of the supersonic annular flow field, boundary con- 
dition and measuring line for physical properties, respec- 
tively. The number of cells is 211x81x31. The adiabatic 
no-slip wall was used as boundary condition. Inflow 
condition was fixed at initial condition. Outflow condi- 


tion is a zero-order extrapolation. Periodic boundary 
conditions were imposed for conservative variables. 
Condensate mass fraction g was set at g = 0 on the wall. 
Table 1 shows calculation conditions used in the 
present calculation. Stagnation pressure po and tempera- 
ture Tọ at stagnation point are 102 kPa and 293 K, re- 
spectively. The initial degree of relative humidity ¢ of 
moist air is 80%. The concentration of the solid particles 
per unit volume net was changed in the range from 1.0 x 
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Fig. 4 Computational domain and boundary conditions 


Table 1 Calculation conditions 
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10° to 1.0 x 10'° 1/m’. The radius of the solid particle R, 
was set to 1.0 x 10° or 1.0 x 10° m and swirl number 
Swin at the inlet was set to 0 or 2.0. 


Results and Discussions 


Figure 5 shows static pressure distributions, which are 
obtained on the basis of the nozzle geometry in Fig.2, on 
the outer wall for dry air (¢@ = 25 %) and moist air (ø = 
70 %) in case of no swirl. In this figure, comparison be- 
tween the experimental and simulated static pressure 
distributions is shown and it is found that simulated re- 
sults agree well with experimental results qualitatively. 

Figure 6 shows contour maps of condensate mass 
fraction g and sonic lines for homogeneous. As seen from 
this figure, positions of sonic line and onset of condensa- 
tion move upstream in case of swirling flow (Sy,in=2.0) 
and this tendency appears particularly in the range close 
to the inner side. 

Contour maps of condensate mass fraction g and sonic 
lines for heterogeneous condensation are shown in Figs. 
7(a) (neto =1.010'° m°, R, = 1.0x10 m) and 7(b) (neto 
=1.0x10"° m”, Ry = 1.0x10° m). In each figure, simu- 
lated results for Sw;in=0 and 2.0 are shown including the 
position of g =1.0x10". It is found from these figures 
that positions of sonic line move upstream in case of 
swirling flow (S\,in=2.0) and this tendency is similar to 
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Fig. 5 Distributions of static pressure on the outer wall 
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Fig.7 Contour maps of g and sonic lines (Mnet, =1.0x 10 m? 


that in case of homogeneous condensation in Fig.6. Fur- 
thermore, positions of the occurrence of heterogeneous 
condensation move upstream largely compared to those 
for homogeneous condensation. Flow characteristics in 


other combination of Meo and R, showed a similar ten- 
dency to that of Fig.7. 

Figure 8 shows distributions of static pressure p/po, 
and nucleation rate Jhom along line A-A’ in Fig.4(c) for Rp 
= 1.0x10° m and 1.0*10° m (neto =1.0*10"° m”, Sw,in=0 
and 2.0). In this figure, results for dry air and homoge- 
neous condensation are also shown for reference. As seen 
from this figure, static pressures for swirling flow be- 
come small compared with those of no swirl and the po- 
sition of the maximum of Jom in case of Sy,in=2.0 is lo- 
cated upstream of that for Swin=0. Further, the positions 
that pressure begins to deviate from an isentrope (dry air) 
for heterogeneous condensation are almost the same as 
one of homogeneous condensation. These distributions in 
other combination of peo and R, showed a similar ten- 
dency to that of Fig.8. 

Effect of the radius of the solid particle R, on distribu- 
tions of condensate mass fraction Znom, Zhet and g (= 
LhomtZhet) along line A-A’ in Fig.4(c) is shown in Fig.9 
for cases of Syin=O and 2.0 (Aneto =1.0x10 m°). It is 
found from this figure that the occurrence position of 
heterogeneous condensation is located upstream com- 
pared with that of homogeneous condensation and g for 
Sw,in=0 or 2.0 is almost the same as one of homogeneous 
condensation regardless of R, at nozzle exit (x/D =3.0). 
Distributions of g in other combination of peo and R, 
showed a similar tendency to that of Fig.9. 
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Fig. 8 Distributions of static pressure and nucleation rate 
(Pinero =1.0*10'° m°) 
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Fig. 9 Distributions of condensate mass fraction (Mnet, = 1.0 x 
10° m°) 


In order to investigate the amount of condensate in the 
range of the outer wall side in the annular nozzle, a new 
function defined by the following equation was intro- 
duced. 


outer Y outer 
G= 2agpurdr / f i 2agpurdr (3) 


vat inner 

where G (radial distribution ratio of liquid phase) 
shows ratio of integrated value of mass flow rate of the 
liquid phase from y=r; (position of the radius corres- 
ponding to half of annular area) to outer wall to inte- 
grated value of mass flow rate of the liquid phase from 
inner wall to outer wall in the cross section in the annular 
nozzle channel perpendicular to the x-axis. 

Distributions of G along x-axis for heterogeneous 
condensation are shown in Figs. 10(a) (neto =1.0x10? 
m”) and 10(b) (neo =1.0x10'° m°). Results for homo- 
geneous condensation are also illustrated for reference. In 
both figures, simulated results for S,;,=0 and 2.0 are 
shown in cases of R, = 1.0x10° m and 1.0x10° m. In 
downstream of the annular nozzle, G in case with swirl 
for each figure becomes large in comparison with cases 
of no swirling flows. This means that the condensate 
gathers toward the outer wall side compared to cases of 
no swirling flows. Further, G for S,,;,=0 or 2.0 is almost 
the same as one of homogeneous condensation at nozzle 
exit (x/D =3.0) regardless of value of Rp. Distributions of 
G in other combination of Meo and R, showed a similar 


tendency to that of Fig.10. 

Figure 11 shows distributions of displacement thick- 
ness of boundary layer 6* on the inner and outer walls 
(Mnet,0 =1.0x10'° m^’). From this figure, it is found that the 
displacement thickness in the case with swirl becomes 
thick in the inner wall side and thin in the outer wall side 
in comparison with the case of no swirl. Further, 6* for 
Syin=0 or 2.0 is almost the same as one of homogeneous 
condensation at nozzle exit (x/D =3.0) regardless of value 
of R, Distributions of 6* in other combination of Mpet,o 
and R, showed a similar tendency to that of Fig.11. 

Figure 12 shows distributions of kinematic viscosity v 
in cross-section at x/D = 3 (reto =1.0x10'° m^’). As seen 
from this figure, v for Sy,i,=0 or 2.0 is almost the same as 
one of homogeneous condensation regardless of R, at 
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Fig. 10 Distributions of G 
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Fig. 11 Distributions of displacement thickness (peto = 1.0 x 
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Fig. 12 Distributions of kinematic viscosity (x/D=3, peto = 
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nozzle exit and the kinematic viscosity becomes large on 
the inner wall side in case with swirl compared with that 
of no swirl. As a result, it is considered that displacement 
thickness becomes thin at the outer wall side and thick at 
the inner wall side as shown in Fig.11. Distributions of v 
in other combination of Meo and R, showed a similar 
tendency to that of Fig. 12. 

As is evident from Figs.8 to 12, the change in the 


physical properties due to heterogeneous condensation is 
similar to that due to homogeneous condensation. As a 
result, it was found in the present calculation condition 
that homogeneous condensation can predict flow proper- 
ties by heterogeneous condensation to some extent. 


Conclusions 


A computational study has been made to investigate 
the effect of the heterogeneous condensation on the 
characteristics of the swirling flow field in a supersonic 
annular nozzle and the differences between homogeneous 
condensation and heterogeneous condensation in the flow 
field. The results obtained are as follows: 

(1) The generating position of condensate in the annu- 
lar nozzle in case of heterogeneous condensation moved 
upstream compared with the case of homogeneous con- 
densation. 

(2) The position that pressure begins to deviate from 
an isentrope (dry air) for heterogeneous condensation 
was almost the same as one of homogeneous condensa- 
tion. 

(3) Condensate mass fraction and radial distribution 
ratio of liquid phase for S,,in= 0 or 2.0 is almost the same 
as one of homogeneous condensation regardless of the 
radius of the solid particle downstream of the nozzle. 

(4) The displacement thickness in the case with swirl 
became thick in the inner wall side and thin in the outer 
wall side in comparison with the case of no swirl and 6* 
for S\,in= 0 or 2.0 is almost the same as one of homoge- 
neous condensation regardless of the radius of the solid 
particle downstream of the nozzle. 

(5) The kinematic viscosity became large on the inner 
wall side in case with swirl compared with that of no 
swirl and kinematic viscosity for Sy,in= 0 or 2.0 was al- 
most the same as one of homogeneous condensation re- 
gardless of the radius of the solid particle downstream of 
the nozzle. 

(6) In the present calculation condition, homogeneous 
condensation could predict flow properties by heteroge- 
neous condensation to some extent. 
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